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ABSTRACT 

We present high time-resolution spectroscopy of the ecHpsing cataclysmic variable 
SDSS J143317.78+101123. 3 obtained with QUCAM2, a high-speed/low-noise electron- 
multiplying CCD camera. iLittlefair et al.l ([200^ measured the mass of the secondary 
star in SDSS J143317. 78-1-101123. 3 using a light-curve fitting technique and obtained a 
value of M2 — 0.060 ± 0.003 Mq, making it one of the three first bona-fide detections 
of a brown-dwarf mass donor in a cataclysmic variable. In this paper we present a 
dynamical measurement supporting this important result. We measured the radial- 
velocity semi-amplitude of the white dwarf from the motion of the wings of the Ha 
emission line and obtained a figure of Ki = 34 ± 4 km s"""^, in excellent agreement 
with the value of if i = 35 ± 2 km s~^ predicted by Littlefair et al.'s model. 
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1 INTRODUCTION 

, Cataclysmic variables (CVs) are close binary stars in which 
• a white-dwarf (WD) primary accretes material from an ac- 
cretion disc f ed by a red-dwarf secondary (see jjollicr (2001) 
, and fWarneii (|l995h for reviews). Theory predicts that the 
■ mass transfer from the secondary to the primary is accom- 
panied by a shortening of the orbital period. This is be- 
lieved to continue until the secondary star drops below the 
hydrogen-burning limit (~ O.O8M0) and becomes a degen- 
erate, brown-dwarf like object, at which point the period 
begins to increase. Models show that up to 70% of all CVs 
in the Galaxy should have brown-dwarf secondary stars at 
the present time, yet not a single such obje ct had been 
positiv ely identified until the discovery by iLittlefair et al.l 
(|2006l ) of a secondary star of mass M2 = 0.052 ± 
O.OO2M0 in the CV SDSS J103533.03-h055158.4. Since then, 
tw o more brown- dwarf m ass donors have been discovered 
by ILittlefair et al. ' (2008) SDSS J150137.22-H550123.3 and 
SDSS J150722. 30-1-523039.8 and it now appears as if the 
missing population of post-period minimum CVs has finally 

been identified. 

The results of lLittlefair et all l|2008l) and lLittlefair et all 
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l|2006l ) rely on an eclipse light-curve fitting technique applied 
to broad-band photometric data. The assumptions underly- 
i ng th is technique appear to be robust (see iLittlefair et al.l 
(j2003) for a discussion), but the results are of such signifi- 
cance that it is important they are independently verified. 
In this paper we present time-resolved spectroscopy of the 
CV SDSS J143317.78+101123.3 (hereafter SDSSJ1433), ob- 
tained with the aim of measuring the radial-velocity semi- 
amplitude of the white dwarf (Ki) and comparing this di- 
rect, dynamical measur ement of Ki wi t h tha t predicted by 
the light-curve model of lLittlefair et af] (|2008l ). The observa- 
tions were obtained with an electron-multiplying CCD (EM- 
CCD), to the best of our knowledge the first time that such 
a device has been used for astronomical spectroscopy. 



2 OBSERVATIONS AND REDUCTION 

SDSSJ1433 is a challenging target as it is faint {g' — 18.5) 
and has a short period (P=78.1 min). This means that spec- 
tra taken with a conventional CCD would be swamped by 
read noise and for this reason we chos e to use an Electron 
Multiplying CCD (|Mackav et al.ll200ll ). These devices are 
able to amplify the signal to such an extent that it renders 
the readout noise negligible. We observed SDSSJ1433 using 
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Table 1 . Full- width at zero intensity, full- width at half maximum, 
peak-to-peak separation and equivalent width of the Ha emission 
line shown in Figure [T] The radial velocity semi-amplitude of the 
WD and the systemic velocity derived in Sections 13. 41 and 13. 5l are 
also listed. 



FWZI 5000 ± 500 km/s 

FWHM 2200 ± 200 km/s 

Peak Separation 1300 ± 200 km/s 

Equivalent Width 147.2 ± 0.5 A 

Ki 34 ± 4 km/s 

7 75 ± 10 km/s 



QUCAMfl an EMCCD-based camera mounted on the red 
arm of the ISIS spectrograpl{3 on the 4.2m WiUiam Herschel 
Telescope, La Palma. QUCAM2 employs a Ik x Ik CCD201 
detector manufactured by E2V which has a frame-transfer 
architecture, providing a dead time of only 12 ms between 
each frame. With the low predicted K\ it was decided to use 
the highest dispersion grating available, the R1200R. This 
provided a wavelength range of 6480-6700A at a spectral 
resolution of 32 km s~^ (0.7A) with a slit width of 1". 

On 2008 April 16, we obtained 628 spectra, each of 30s 
exposure time. The moon was full, the seeing was approx- 
imately 0.7-1" and the sky transparency was variable. A 
nearby comparison star was also placed on the slit to correct 
for slit losses and changes in transparency. Arc spectra were 
taken every 40 minutes for wavelength calibration. No flux 
standard was observed. This is justified because there was 
intermittent cloud, we only observed over a narrow wave- 
length range and we were only interested in extracting ve- 
locities from the data. 

The spectra were extracted using a simple spatial bin 
across three windows within each frame covering the target, 
the reference star and the large area of sky in between. The 
sky-subtracted s pectra were cas t into 40 phase bins using 
the ephemeris of iLittlefair et ahl l|2008l ). and then corrected 
for slit losses and transparency variations by dividing by the 
integrated fiux in the comparison star spectra. 



3 RESULTS 

3.1 Averaged spectrum 

Figure [1] shows the average spectrum of SDSSJ1433. The 
most prominent feature is the broad, double-peaked emis- 
sion line of Ha. The equivalent and velocity widths of the 
line are given in Table[T] A weaker feature due to He I A6678 
is also visible. The average spectrum is typical of other 
eclipsing dwarf novae below the period gap, e.g. WZ Sge 
jSkidmore et al.ll2000l ). 

3.2 Continuum and Ha light curves 

We computed the light curve of the continuum by summing 
the flux in line-free portions of each spectrum. We then fit- 
ted and subtracted the continuum of each spectrum, and 

^ http:/ /www. ing.iac.es/Engineering/detectors/qucam2. html 
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Figure 1. The average spectrum of SDSSJ1433, obtained by com- 
bining all 628 spectra. 
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Figure 2. Light curves of the continuum only and the Ha emis- 
sion line. 

summed the residual flux in the Ha emission line. The re- 
sulting light curves are shown in Figure [2] 

Figure [2] shows that Ha experiences broader eclipses 
than the continuum. The eclipse depth for Ha is ~ 50 %, 
while for the continuum it is ~ 75 %. This is consistent with 
Ha being emitted in the optically thin outer parts of the 
accretion disc and the hotter, optically thick, inner parts of 
the disc being responsible for the majority of the continuum 
emission. This picture is confirmed by analysing the light 
curves of the wings and core of the line separately, which 
show a deeper eclipse in the former than the latter. An or- 
bital hump prior to eclipse is also apparent in the Ha light 
curve of figure [2] caused by the changing aspect of emission 
from the bright spot. 

3.3 Trailed spectra and Doppler tomogram 

The phase-binned, continuum-subtracted Ha profile is 
shown as a trail in the lower panel of Figure |3] The primary 
eclipse around phase is clearly seen, as is the rotational dis- 
turbance where the blue peak of the emission line is eclipsed 
prior to the red peak. The bright spot is clearly visible as an 
S-wave moving between the two peaks. The orbital modula- 
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tion of the wings of the emission hne can arguably just be 
made out. There is also evidence for a 'shadow' on the blue 
edge of the S-wave around phase 0.25 which appears remark- 
ab ly similar to the models of gas stream overflow computed 
bv lHellier fc RobinsonI l|l994) . 

We computed the Doppler tomogra m of th e Hq trail 
using Fourier-filtered back projection (see iMarshI l|200lh for 
a review). As expected, the dominant feature is the ring of 
emission representing the accretion disc which appears to 
be centred on the expected position of the WD. The bright 
spot lies on this ring at a velocity intermediate between the 
free-fall velocity of the gas stream and the Keplerian disc 
velocity along its path, as o bserved in som e other CVs, e.g. 
U Ge m l|Marsh et al.lll990l ) and WZ Sge (jSpruit fc RuttenI 
Il998h . There is no evidence for emission from the secondary 
star, nor for asymmetries in the disc emission due to, for 
example, spiral shocks. 

3.4 The diagnostic diagram 

The WD primary star orbits with velocity Ki around the 
centre of mass of the system, resulting in a periodic Doppler 
shift of any light it emits. If we imagine the simplified case 
of an ideal accretion disc centred on the WD with a per- 
fectly symmetric brightness distribution, then the velocity 
shift of the emission lines generated in the disc will be sub- 
ject to the same modulation as the WD. If we measure the 
line centroids and plot them as a function of orbital phase 
then we should see a sine wave corresponding to the orbital 
motion of the WD. More realistically, however, if we include 
the light emitted by the bright spot, we can see that calcu- 
lated line centroids will, rather than precisely following the 
WD, be perturbed in the direction of the bright spot. The 
resultant radial velocity (RV) curve will not only have an 
excessive amplitude but will also be offset in phase with re- 
spect to the true motion of the WD. It is thus important to 
exclude the bright spot contribution by only examining the 
light emitted in the wings of the line profile. These emissions 
correspond to material orbiting within the disc at small radii 
from the primary star, where they are only minimally con- 
taminated by the bright spot. If the RV curve that we obtain 
using this technique has the correct phase offset relative to 
the primary eclipse then we can be confident of an accurate 
semi-amplitude measurement. 

To measure the radial velocities of the wings of the 
H q emission line we u sed t he double-Gaussian technique 
of ISchneider fc Youngj (|l98(]| ). The best signal-to-noise ra- 
tio was obtained using a Gaussian width of 350 km s^^ 
and the separation of the Gaussians was then varied be- 
tween 1500 and 3000 km s~^ in steps of 100 km s~^ so as 
to explore annular regions of the accretion disc of decreas- 
ing radii. At each Gaussian separation the resulting radial 
velocities were fitted with a sine function. The phase offset 
((^o), semi-amplitude (K), systemic velocity (7) and frac- 
tional error in the amplitude {ok/K) of the sine fits were 
then plotted against the Gaussi an separation as a diagnos- 
tic diagram (IShafter et al.l[l986t ). as shown in Figure |4] An 
example of the RV data and sine fit is shown in Figure [5l 
for a Gaussian separation of 2900 km s~^. 

The standard way of obtaining Ki from a diagnostic di- 
agram is to use the value corresponding to where the / K 
curve is at a minimum. If we do this we find a value of 
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Figure 3. Bottom: the Ha trailed spectrum, with two cycles plot- 
ted for clarity. Top: the Doppler map of Ha. The three crosses rep- 
resent the centre of mass of the secondary star (upper cross), the 
system (middle cross) and the WD (lower cross). The Roche lobe 
of the secondary star, the predicted gas stream trajectory and the 
Keplerian velocity of the disc along the stream are also shown, cal- 
culated assuming the system parameters given in iLittlefair et aP 
12008). The circular tick marks represent steps of 0.1 Li (the 
inner Lagrangian distance) towards the WD. The asterisks rep- 
resent the points of closest approach to the WD. 



Ki = 44km s~^. However, it can be seen that the corre- 
sponding value of 00 is non-zero and declines towar ds zero 
at hig her Gaussian separations. As pointed out by iMarshI 
(198«), this highlights the main drawback of diagnostic dia- 
grams - the derived iCi-value is noise dependent. The best 
solution is to use those data points corresponding to the 
highest Gaussian separations to extrapolate to a Jfi-value 
whose corresponding 00 = 0. This require s the construction 
of a 'Light Centre' diagram (|Marshlll988l '). 
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Figure 4. The diagnostic diagram for H« - see text for details. 



3.5 The light centre method 

At no point on the diagnostic diagram does the sine fit 
have a zero phase offset relative to the WD, although at 
the extremes of the H a wing s it get s close. The light centre 
method, described bv lMarsh (|l988r ). offers a way of extrap- 
olating the data to estimate what the amplitude of the RV 
curve would be at 00 = 0. Just as the Doppler map is a view 
of an emission line transformed into velocity space, so the 
light centre diagram, shown in Figure [S] is a velocity space 
projection of the RV data plotted in the diagnostic diagram. 



The radial velocity of the WD can be determined from a 
light centre diagram by extrapolating a linear fit to the RV 
data points and reading off the y-axis intercept. To avoid 
contamination from the bright spot we only included the 
nine largest Gaussian separations in the linear fit. The re- 
sulting intercept is Ki = 34 ± 4 km s~^, the error quoted 
corresponding to a combination of statistical error and sys- 
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Figure 5. The measured RV data for a Gaussian separation of 
2900 km s^^ and the resulting sine curve fit (solid line). Only 
data points marked by a square were included in the fit, i.e. the 
rotational disturbance was not taken into account. The dashed 
line shows the sine curve corresponding to the motion of the WD 
as derived from the light centre diagram presented in Figure [6] 
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Figure 6. The light centre diagram for H«, where — i^sin^o is 
plotted on the abscissa and —KcosijiQ on the ordinate. A linear 
fit to the points marked by squares is also shown, extrapolated to 
where it intercepts the y-axis. The smallest Gaussian separation 
corresponds to the top-most point and the largest to the right- 
most point. 



tematic error due to the exact choice of which points to 
exclude from the linear fit. The RV curve corresponding to 
this value of Ki and 00 = is shown as the dashed curve in 
Figure [5] 

We are confident that our measurement of K\ accu- 
rately reflects the motion of the WD in SDSSJ1433 for the 
following reasons. First, the bottom panel of figure |4] clearly 
shows a trend to zero phase offset at large Gaussian separa- 
tions, as one would expect if one is measuring the real radial 
velocity of the WD in combination with a steadily decreas- 
ing contribution from the bright spot. In such instances, the 
light centre techn ique has bee n prove n to give reliable re- 
sults (e.g. IP Peg, lMarshlll988l : U Sco, iThoroughgood et al.l 
l200j ). Second, the top panel of figure |4] shows that the ra- 
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dial velocity semi-amplitude is remarkably independent of This paper has been typeset from a Tp^X/ ffl^jX file prepared 
the Gaussian separation, indicating that our systematic er- by the author, 
rors are indeed small. 



4 CONCLUSIONS 

We have measured the radial velocity amplitude, Ki = 
34±4 km s-\ of the WD in the short-period CV SDSSJ1433 
by analysing the RV motion of the extreme wings of the Hq 
emission line. In order to account for the velocity contami- 
nation from the conspicuous emission of the bright spot, the 
light centre technique was used. The measured value of Ki 
is consistent with the predicted value of K i = 35 ± 2 km 
s~^ fr om the light-curve fitting technique of iLittlefair et aP 
l|2008h . Our result therefore supports the validity and ac- 
curacy of the purely photometric technique of measuring 
masses and argues in favour of the presence of a brown-dwarf 
donor in SDSSJ1433, and by implication in SDSSJ1501 and 
SDSSJ1507 also. The results presented in this paper also 
demonstrate that the combination of a large-aperture tele- 
scope, an intermediate-resolution spectrograph and an EM- 
CCD is uniquely capable of tackling this type of photon- 
starved observation. 
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